Particle motion in a deformed potential using a transformed oscillator basis 
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The quantum description of a particle moving in a deformed potential is investigated. A pseu- 
dostate (PS) basis is used to represent the states of the composite system. This PS basis is obtained 
by diagonalizing the system Hamiltonian in a family of square integrable functions. In this work the 
Transformed Harmonic Oscillator (THO) functions, obtained from the solutions of the Harmonic 
Oscillator using a Local Scale Transformation (LST), are used. The proposed method is applied to 
the n Be nucleus, treated in a two-body model ( 10 Be + n). Both structure and reaction observables 
have been studied. 

Wavefunctions and energies obtained for the bound states and some low-lying resonances are 
compared with those obtained by direct integration of the Schrodinger equation. The dipole and 
quadrupole electric transition probabilities for the low-energy continuum have been calculated in 
the THO basis, and compared with the exact distributions obtained with the scattering states. 
Finally, the method is applied to describe the 11 Be states in the Coulomb breakup of 11 Be+ 208 Pb 
at 69 MeV/nucleon. The energy and angular distributions of the exclusive breakup have been 
calculated using the Equivalent Photon Method, including both El and E2 contributions. The 
calculated distributions are found to be in good agreement with the available experimental data 
from RIKEN [Phys. Rev. C70, 054606]. At the very forward angles, the cross section is completely 
dominated by the dipole couplings. 

PACS numbers: 24.10.Eq, 25.10.+S, 25.45.De 



INTRODUCTION 



It is well known that the quantum collision of a weakly 
bound system by a target is influenced by the coupling 
to the unbound states of the projectile. For nuclear col- 
lisions, this effect was first noticed in deuteron-induced 
reactions, and later observed in the scattering of other 
loosely bound nuclei, such as halo nuclei. Several re- 
action frameworks have been envisaged to account for 
this effect. Among them, one of the most successful 
has been the Continuum-Discretized Coupled-Channels 
(CDCC) method [H, 0, originally developed to account 
for the breakup channels in deuteron scattering and later 
extended to other weakly-bound nuclei, such as 6,7 Li, 
11 Be, or 8 B, among others. In all these cases, the pro- 
jectile is described in a two-body model ( 6 Li= 4 He+d, 
7 Li= 4 He+ 3 H, n Be= 10 Bc+n, etc) and the method con- 
siders explicitly the possible dissociation of the projec- 
tile into its two fragments. In its standard formulation, 
the excitation of each fragment is nevertheless ignored. 
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This is a good approximation for deuteron scattering, 
for which both constituents can be considered inert at 
the energies of interest in nuclear studies, but it is more 
questionable for more complex systems. Moreover, the 
bound and unbound states of the two-body system are 
considered to be well described by pure single-particle 
configurations. This approximation ignores possible ad- 
mixtures of different core states in the wave functions of 
the projectile. These admixtures are known to be impor- 
tant, particularly in the case of well-deformed cores, as 
for example in the n Be halo nucleus. 

A recent attempt to accommodate these effects within 
the CDCC method was done in Ref. Q, and applied 
to the scattering of one-neutron halo nuclei. In that 
work, the states of the corc+valcnce system were de- 
scribed within the particle- rotor model [4j . The unbound 
states of the compound system were described by contin- 
uum bins which, following the standard procedure, were 
constructed by superposition of scattering states. These 
scattering states are obtained by direct integration of 
the Schrodinger equation, with the appropriate bound- 
ary conditions. 

Alternatively, the bound and unbound states of the 
system can be obtained by diagonalizing the Hamilto- 
nian in a suitable basis of square-integrable functions. 
The eigenfunctions of the system are expressed as an ex- 
pansion in the basis functions. In practical calculations, 



the basis needs to be truncated, leading to a finite ex- 
pansion of the eigenfunctions. Therefore, these states 
and their corresponding eigenvalues can be regarded as a 
finite approximation to the exact states of the system and 
are referred hereafter as pseudo-states (PS). This proce- 
dure has been applied, for example, to describe the states 
of two-body nuclei interacting via a central potential [a- 
7] and, more recently, also for three-body nuclei (8l— flljj - 
A variety of bases have been used in these applications, 
such as harmonic oscillator (HO), Gaussian, Laguerre 
functions, etc. The procedure can be also extended to 
deformed systems. A natural choice for the PS would be 
the deformed HO potential |12|,[l3|]. However, this basis 
is not suitable to describe the bound states of weakly- 
bound nuclei due to its Gaussian asymptotic behavior. 
Several alternatives have been proposed in the literature, 
for example, the eigenstates of a truncated Woods-Saxon 
potential [1J| or the Sturmian basis pjl [l6J . 

In this work, we propose the use of a Transformed Har- 
monic Oscillator (THO) basis to describe the states of a 
two-body system mutually interacting with a deformed 
potential. This basis has been previously applied to the 
case of spherical systems Q so we present here its exten- 
sion to deformed systems. The THO basis is obtained 
by applying a Local Scale Transformation (LST) to the 
Harmonic Oscillator (HO) basis. The LST. adopted from 
a previous work of Karataglidis et al. [17[, is such that 
it transforms the Gaussian asymptotic behavior into an 
exponential form, thus ensuring the correct asymptotic 
behavior for the bound wave functions. The accuracy 
of this THO basis was tested for several reactions in- 
duced by deuteron and halo nuclei, showing an excellent 
agreement with the standard binning method, and an im- 
proved convergence rate in the case of narrow resonances 

to 

For a deformed potential, the calculation of bound 
and unbound states becomes a multi-channel problem, 
since, in general, for each physical state there will be 
contributions from several orbital angular momenta and 
core states. For bound states, the calculation of the 
energies and eigenfunctions is analogous to the single- 
channel case, because these quantities are directly ob- 
tained from the diagonalization in the chosen PS basis. 
For unbound states, the eigenfunctions (and their cor- 
responding eigenvalues) obtained from the Hamiltonian 
diagonalization can be regarded as a finite and discrete 
representation of the exact states. In general, resonances 
(quasi-stationary states) correspond to combinations of 
these positive-energy eigenstates and hence its identifi- 
cation is not straightforward. 

For a particle moving in a central potential (with pos- 
sibly a spin-orbit component) this is a relatively straight- 
forward problem and indeed a variety of methods have 
been proposed to compute resonance energies and widths. 
For example, they can be obtained from the poles of the 
S'-matrix in the complex energy plane. A simpler method 
is to define the resonance as the energy at which the 
phase-shift crosses tt/2. The width is then obtained from 



the inverse of the derivative of the phase-shift, evaluated 
at the energy of the resonance. These methods rely on 
the knowledge of the scattering states at large distances 
(from which the S'-matrix and hence the phase-shifts can 
be extracted) and then they cannot be directly applied 
to PS methods, given the wrong asymptotic behaviour 
of the PS functions. In this case, the identification of 
resonances can be done using the so-called stabilization 
method [19|, [2CJ • This is a procedure envisaged to iden- 
tify and construct the most localized continuum wave 
functions when the positive energy states are expanded 
in a discrete basis, depending on one or more parame- 
ters. In practice, this can be achieved by diagonalizing 
the Hamiltonian as a function of these parameters (for 
example, the basis size) and then scanning the resultant 
eigenvalues for the continual appearance of a stabilized 
value which, unlike the others, is insensitive to the size 
of the basis. In some previous works, we have success- 
fully applied this technique to obtain the resonances of 
two-body systems with central potentials using the THO 
basis 0, [l8[- In this work, we explore the validity of 
this method for the multi-channel situation that arises in 
the deformed case. Our aim with this work is to assess 
the capability of the THO basis for calculations includ- 
ing core deformation in the simpler two-body systems, as 
n Be. This step is necessary and unavoidable for provid- 
ing a solid basement to proceed with the generalization 
of the formalism to more challenging situations, such as 
the case of three-body composite systems including core 
deformation, or to the scattering of a two-body system 
by a third body, including core deformation in one of the 
clusters of the composite system. 

The work is structured as follows. In Section [XT] the 
THO method based on the parametric LST is reviewed 
and the structure model used in subsequent calculations 
is discussed. In Section IIII1 general expressions for the 
electric transition operators for the particular case of a 
two-body system with a deformed core are provided. In 
Sec. IIVI the model is applied to describe the structure 
of the 11 Be nucleus. The basis so obtained is then used 
to describe the Coulomb breakup of 11 Be on 208 Pb at 
69 MeV/nucleon, comparing our results with the avail- 
able data. Finally, in Section [V] the main results of this 
work are summarized. 



II. EIGENSTATES OF A DEFORMED 
POTENTIAL IN A PS BASIS: THE THO BASIS 

In this section, we briefly review the features of the 
PS basis used in this work. This basis is an exten- 
sion of the THO basis used in our previous works to 
describe the states of a composite system consisting of 
two interacting inert fragments, such as a valence par- 
ticle (proton/neutron) and a spherical and stable core. 
The goal of this extension is to allow core-excited admix- 
tures in the description of the states of the composite 
system and hence the possibility of dynamic core excita- 



tion mechanisms in reactions involving these nuclei. For 
completeness, we review first the situation in which the 
core degrees of freedom are neglected. In this case, the 
core+valence Hamiltonian is simply given by: 

H = T r + V vc (r) (1) 

where r is the relative coordinate between the valence 
and the core, T r the core- valence kinetic energy operator 
and V vc (r) is the interaction between the valence particle 
and the core. The eigenstates of this Hamiltonian can 
be characterized by the energy eigenvalues (e) and the 
set of quantum numbers {£,s,j}, which correspond to 
the orbital angular momentum (£), the valence spin (s) 
and their sum (j = £ + £)• For a central potential with, 
possibly, a spin-orbit term, these states can be written 
as: 

4>e,t,jif) = Re,l,j{r)yisjm{f) (2) 

where yg sjm (f) = [Ye(r) <g> Xs] jm , with Xs a spin func- 
tion. The radial functions R e> ij[r) can be obtained by 
solving the Schrodinger equation subject to the appro- 
priate boundary condition for bound (s < 0) or un- 
bound (e > 0) states. Alternatively, these functions 
can be obtained by diagonalizing the Hamiltonian ([T]) 
in a discrete basis. Since any complete basis will be in- 
finite, this procedure is not feasible in practice unless 
the basis is truncated. By doing so, one obtains a fi- 
nite (an approximate) expansion of the functions R(r) in 
the selected basis. If the basis functions are denoted by 
¥>n,i,j(r) = XnA r )yesj7n(r), we will have: 



N 



R p{ r ) = ^Cfi,nXnA r ) 



(3) 



n=l 



where j3 = {e,£,s,j} and N is the number of states re- 
tained in the basis. 

As already mentioned, there are many possible choices 
for the basis functions {if n } (Gaussians, harmonic oscil- 
lator, Laguerre, etc). In this work we use the transformed 
harmonic oscillator (THO) basis, obtained from the har- 
monic oscillator basis with an appropriate LST [2l|, [22| ■ 
If the LST function is denoted by s(r), the THO states 
are obtained as 



R. 



THO 




(4) 



where R^f{s) is the radial part of the HO functions. 
With the criterion given above, the LST is indeed not 
unique. In Ref. [23| the LST was defined in such a way 
that the first HO state is exactly transformed into the 
exact ground state wave function, assuming that this is 
known. Therefore, by construction, this wave function is 
exactly recovered for any arbitrary size of the basis. In a 
more recent work Ml we adopted the parametric form of 
Karataglidis et al. [13] 



s(r) 



V2b 



\r) + ( 7 V? 



(5) 



that depends on the parameters m, 7 and the oscilla- 
tor length b. Note that, asymptotically, the function 
s(r) behaves as s(r) ~ jfv^3 anc ^ hence the functions 
obtained by applying this LST to the HO basis behave 
at large distances as exp(— r y 2 r/2b 2 ). Therefore, the ra- 
tio 7/6 can be related to an effective linear momentum, 
k c a = 7 2 /26 2 , which governs the asymptotic behaviour of 
the THO functions. As the ratio 7/6 increases, the radial 
extension of the basis decreases and, consequently, the 
eigenvalues obtained upon diagonalization of the Hamil- 
tonian in the THO basis tend to concentrate at higher 
excitation energies. Therefore, 7/6 determines the den- 
sity of eigenstates as a function of the excitation energy. 
In all the calculations presented in this work, the power 
m has been taken as m = 4. This choice is discussed 
in Ref. [17] where the authors found that the results are 
weakly dependent on m. 

Note that, by construction, the family of functions 
Rnt( r ) are orthogonal and constitute a complete set 
with the following normalization: 



r a |i£r(r)| a dr = l. 



(6) 



Moreover, they decay exponentially at large distances, 
thus ensuring the correct asymptotic behaviour for the 
bound wave functions. In practical calculations a finite 
set of functions (j4|) is retained, and the internal Hamil- 
tonian of the projectile is diagonalized in this truncated 
basis with N states, giving rise to a set of eigenvalues and 
their associated eigenfunctions, denoted respectively by 
{e„} and {</r e (r)} (n = 1, . . . , N). As the basis size is 
increased, the eigenstates with negative energy will tend 
to the exact bound states of the system, while those with 
positive eigenvalues can be regarded as a finite represen- 
tation of the unbound states. 

The formalism can be extended to the situation in 
which the core degrees of freedom are taken into account 
explicitly In this case, the Hamiltonian (TTJ) is generalized 
to 



H = T r + V vc {f,0 + h COTS (0 



(7) 



where /i C ore(£) is the intrinsic Hamiltonian of the core, 
whose eigenstates will be denoted by {<j>iMi}- Addi- 
tional quantum numbers, required to fully specify the 
core states, are not included for notation simplicity. Note 
that the valence-target interaction, V vc (r, £), contains 
now a dependence on the core degrees of freedom (de- 
noted generically by £). 

The eigenstates of the Hamiltonian cannot any longer 
be written in the form of Eq. @. Instead, these states 
will be a superposition of several valence configurations 
and core states, i.e. 



* 



e:JM 



(?,c)=Eu) 



yisjm(r) (8) <j>i(£) 



JM 



(8) 



Upon replacement of the expansion (jSJ into the 
Schrodinger equation, one gets a coupled set of differ- 
ential equations for the radial functions R ea (r). For 



bound states, these radial functions decay exponentially 
for r — > oo giving rise to square-integrable functions. For 
continuum states, the functions R E , a (r) are also obtained 
by solving a set of coupled radial equations, but subject 
to the boundary condition that incident waves occur only 
in the entrance channel characterized by a given set of 
quantum numbers a = {£, s, j, I}. Therefore, for each 
continuum energy, there are as many scattering solutions 
as possible values of a, compatible with the total angular 
momentum J. 

Alternatively, the functions R £ , a (r) can be obtained 
using an expansion in a PS basis, such as the THO basis 
described above. In this case, the basis must include also 
the new core degree of freedom 



K,jm(t,0 = K H a u {r) y esjm (f) ® MO 



JM 



(9) 



In this basis, the states of the system will be expressed 



as 



iV 



^%(r,i) = EE c «/n,«(^), (10) 



which, for simplicity, is taken to be axially symmetric. 
Thus, we can characterize the deformation by a single 
parameter 02- In the body-fixed frame, the surface ra- 
dius is then parameterized as -R(£) = -Ro[l + /32^2o(0]j 
with Rq an average radius. Starting from a central poten- 
tial, Vvc (r) , the full valence-core interaction is obtained 



by deforming this interaction as, 

V vc (r,0 = V v ^(r-S 2 Y2o(i)), 



(13) 



with §2 = P2R0 being the deformation length. Trans- 
forming to the space-fixed frame of reference, and ex- 
panding in spherical harmonics, this deformed potential 
reads 

V vc {r, () = E V ™ (r)YcM {?)Y£ M (|) , (14) 



C,M 



AC), 



where the radial form factors V vc (r) are obtained by 
projecting the deformed potential (fT3f onto the required 
multipoles. 



where i is an index that labels the order of the eigenstate. 
These eigenstates are spread in the energy spectrum 
with a density strongly related to the basis parameters, 
mainly N and 7/6, and to the continuum structure for the 
selected Hamiltonian, i.e. presence of resonances or dif- 
ferent breakup thresholds. Moreover, this density reflects 
the momentum distribution of the eigenstates which be- 
comes important to obtain continuous energy or momen- 
tum distributions of different observables from their dis- 
crete representation in the PS basis 0, 0, Ojl, [24| . Gener- 
alizing the expression in [l8| , the density of states is here 
defined as: 



N n a 



p(*) = EE^ J /i*ii)' 



(ii) 



where \k a Jf) denotes the exact scattering wavefunction 
for an incoming wave in the a channel. Note that the dif- 
ference between k and k a relays on the threshold energy 
for each channel. 

With this definition the integral of the density with re- 
spect to the momentum is the number of THO functions 
selected (N) times the number of channels (n a ): 



p(k) dk = Nn a , 



(12) 



assuming that we have included N THO functions for 
each channel a. Note that this integrated density is in- 
dependent of the LST parameters. 

The afore-mentioned method can be applied to any 
Hamiltonian of the form ([7]). In the calculations pre- 
sented in this work, the composite system is treated 
within the particle-rotor model 4]. Therefore, we as- 
sume that the core nucleus has a permanent deformation 



III. ELECTRIC TRANSITION PROBABILITIES 
IN THE PS BASIS 

The accuracy of the PS basis to represent the con- 
tinuum can be studied by comparing the ground-state to 
continuum transition probability due to a given operator. 
Here we consider the important case of the electric dis- 
sociation of the initial nucleus into the fragments c + v. 
This involves a matrix element between a bound state 
(typically the ground state) and the continuum states. 

The electric transition probability between two bound 
states \Ji) and \Jf) (assumed here to be unit normal- 
ized) is given by the reduced matrix element (according 
to Brink and Satchler convention [251 ]) 



B(EX;i^ f) = H±i \ { j f \\M(EX)U)\ 2 . 



(15) 



where A4 is the multipole operator. In a core+valence 
model, the electric transition operator can be written as 
a sum of three terms [26[ : one for the excitation of the 
valence particle outside the core, one for the excitation of 
the core as a whole and one for mixed excitations involv- 
ing simultaneous excitations of core and valence particle, 

A-l k 

M(EX^) = J2 E h( k > m >v) 

fc— 1 m— — k 

x M sp {Ekm)M C ore{E(\- k)(n — mj) 

+ MspiEX^+McoreiEXfJ,), (16) 

where f\(k, m, (/,) is a well-defined function of its indices 
and the single particle contribution has the usual form, 



M sp (EXfi) = Z$er x Y Xli {r), 



(17) 



with the effective charge: 



IV. TEST EXAMPLE: APPLICATION TO U BE 



/(A) 



z,, 



rn v + m c 



Z c 



(18) 



In the case of a transition to a continuum of states, 
\kJf), the definition (IT51) is replaced by (see for example 

my- 

,mEx^2j f + i , vc t KkJfimEX)m ^ (19) 



ds 



2J t + 1 (2ir) 3 h 2 



with fc = ^/2fib c e/h. Note that the extra factor appearing 
in Eq. (|19l) with respect to Eq. (fl5l) is consistent with 
the convention (kJ\k'J) = S(k — k 1 ) and the asymptotic 
behaviour, 



u a i yk a i , rj 






o a i a il[ (k a r) 



: (J) 



~T SX a H v (ka>r) , (20) 



where u a (k a ,r) — R a (k a ,r)r (using an obvious notation 
where the continuum e label has been replaced by a de- 
pendence on the corresponding momentum k). 

Using a finite basis, one may calculate only discrete val- 
ues for the transition probability. According to Eq. (|T5|) , 
the B(E\) between the ground state (with angular mo- 
mentum Ji) and the n-th PS is given by 



B {N) {E\;g.s. -m) = 



2J/ + 1 

2J t + l 



^ n N U\M(EX)\\%. 



g.s. 



(21) 

In order to relate this discrete representation to the con- 
tinuous distribution (| 19 j) one may derive a continuous 
approximation to (|19j) by introducing the identity in the 
truncated PS basis, i.e. 



A' 



AN) 
l JM 



I* 



n,JM/\^n,.JM\ 



(22) 



For N — » co this expression tends to the esaci identity 
operator for the Hilbert space spanned by the eigenfunc- 
tions of the considered Hamiltonian. By inserting (|22|) 
into the exact expression (|19[) we obtain the approximate 
continuous distribution, 

dB(EX) 2J f + 1 n vc k 



1 (2tt) 3 ^ 2 
Y,( kJ f\*n N i)(*l%\\M(E\)\\% 



2J, 

N 



(23) 



This approach provides a smoothing procedure to 
extract continuous distributions, as a function of the 
asymptotic energy e (or, equivalently, the linear momen- 
tum k) , from the discrete distributions obtained with the 
PS basis [3, |28|. This is particularly convenient in situa- 
tions in which the calculation with the scattering states 
themselves is not possible, such as in the CDCC method. 



A. Energy spectrum and wave functions in the PS 

basis 



As an illustration of the formalism presented in the pre- 
ceding section, we consider the 11 Be nucleus. This choice 
is motivated by the fact that this nucleus is one of the 
best known one-neutron halo nuclei. Many of its prop- 
erties can be understood in a simple two-body model, 
comprising a valence neutron orbiting a 10 Be core. For 
example, the ground state (l/2 + ) and the only bound ex- 
cited state (1/2 - ) are reasonably well described by 2s\/2 
and lpi/2 single-particle configurations, relative to the 
10 Be(g.s.) core. Excited states in the continuum are also 
reasonably well described in terms of single-particle ex- 
citations of the halo neutron outside the 10 Be(g.s.) core. 
This single-particle picture has been extensively used in 
the literature to explain also reactions induced by this 
nucleus (see for instance [IS, [29|, l30|). However, there 
are also numerous experimental and theoretical evidences 
that these low- lying states of n Be contain significant 
admixtures of core-excited components [31H341 ] . Conse- 
quently, an accurate description of reactions involving 
this nucleus requires the inclusion of its states beyond 
the simple single-particle picture. 

In the calculations presented in this work, we use the 
particle- rotor model of Bohr and Mottelson with the n Be 
Hamiltonian of Ref. [l6| (model Bel2-b), which consists 
of a Woods-Saxon central part, with a fixed geometry 
(R = 2.483 fm, a = 0.65 fm) and a parity-dependent 
strength. The potential contains also a spin-orbit part, 
whose radial dependence is given by the derivative of the 
same Woods-Saxon shape, and strength V so — 8.5 MeV. 
For the 10 Be core, this model assumes a permanent 
quadrupole deformation /32=0. 67. Only the ground state 
(0+) and the first excited state (2+, E x = 3.368 MeV) 
are included in the model space. For the valence-core 
orbital angular momentum, we consider the values (. < 3. 

To generate the THO basis we use the LST of Eq. © 
with m = 4, b = 1.6 fm and 7 = 1.84 fm 1 / 2 . The value of 
b was determined in order to minimize the ground state 
energy of n Be in a small THO basis. The factor 7/6 leads 
to a k e ff compatible with a maximum excitation energy 
of about 10 MeV, which is enough for the calculations 
presented below. 

Once these parameters have been fixed, the THO ba- 
sis is generated for different values of N, the number of 
oscillator functions, and the convergence of different ob- 
servables is studied with respect to this number. We 
should remark that the total number of basis functions is 
this number times the number of channels n a . However, 
the latter depends on the total angular momentum J of 
the state under consideration, and will be the same in 
any method based on the angular momentum expansion 
of the wave functions. Therefore, we will refer to N as the 
basis size as we understand it is the most honest way of 
comparing with other methods. We find that the ground- 
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FIG. 1: (Colour online) Radial parts of the ground state wave 
function for the n Be nucleus obtained by direct integration of 
the Schrodinger equation (solid lines) and by diagonalization 
in a THO basis with N=15 states (dashed lines). 




FIG. 2: (Colour online) Radial parts of the l/2~ excited 
bound state wave function for the Be nucleus obtained by 
direct integration of the Schrodinger equation (solid lines) and 
by diagonalization in a THO basis with N — 15 states (dashed 
lines) . 



state energy is already fully converged with a relatively 
small basis (N « 15). 

Within the model space used in our calculations (I = 
0,2, £ < 3), there are n a — 3 channels contributing to 
the ground state wave function, namely | 10 Be(0 + )(8>Si 



1 10 



'1/2/ 



Be(2 + )(g)d 3/2 ) and | 10 Be(2+)<g)d5 /2 ). In Fig. [Q we de- 
pict these radial parts of the ground-state wave function 
obtained from the diagonalization of the Hamiltonian in a 
THO basis with N=15 oscillator functions (dashed lines). 
For comparison, we include also the solutions obtained 
by direct integration of the Schrodinger equation (solid 
lines). Both calculations give basically identical results. 
It can be seen, as expected, that the | 10 Bc(0 + ) (g) s 1 / 2 ) 
component is the dominant one, accounting for about 
80% of the norm. This radial component exhibits a node, 
due to the presence of a Pauli forbidden state (arising 
from the lsi/2 orbital in the spherical basis). 

The assumed Hamiltonian reproduces also the position 
of the bound excited state at E x = 320 kcV (l/2~). In- 
deed, this state appears also in the diagonalization of the 
THO basis. The separation energy is reproduced within 
a few percent with a basis of N=15 states and the radial 
components are also found to be in perfect agreement 
with those obtained by direct integration of the coupled 
differential equations. This is shown in Fig. [5] 

We proceed to discuss now the description of reso- 
nances in the PS basis. As explained in the introduction, 
the identification of the resonances is done using the sta- 
bilization method of Hazi and Taylor [ljj, [2CJ , extended 
to the multi-channel case. The procedure is the same as 
in the single-channel case, i.e., we diagonalize the Hamil- 
tonian over either a successively larger basis set or as 
a function of a continuous parameter which defines the 
basis for a given N value. Then, the evolution of the 
spectrum as a function of N or the continuous param- 
eter is studied. When a resonance is present, there are 
some eigenvalues whose energies are stabilized for a range 



of values of N or the continuous parameter. This prop- 
erty has been employed empirically in many works, and a 
formal justification has also been provided by Lippmann 
and O'Malley J3|. 

The selected Hamiltonian contains low-lying reso- 
nances ate = 1.2 MeV (5/2+), 2.7 MeV (3/2-) and 
3.2 MeV (3/2 + ) [lg|. These values are confirmed apply- 
ing the stabilization method with the THO basis, in the 
two ways described above. As an example, in Fig. [31 we 
show the results for J T = 5/2+ . In the upper panel, the 
sequence of continuum states with J 71 ' = 5/2+ is plot- 
ted versus the continuum parameter 7 of the LST, and 
for a fixed value of N (N=10). In the lower panel, the 
J v = 5/2+ eigenvalues obtained from the diagonaliza- 
tion of the assumed Hamiltonian in the THO basis are 
plotted as a function of the discrete basis size parameter 
(N), with 7 fixed to 1.84 fm 1 / 2 . The dashed line marks 
the known location of the first 5/2 + resonance deduced 
from the behavior of the phase-shifts and the dotted line 
marks the n+ 10 Be(2 + ) threshold. In both plots, the en- 
ergy stabilization precisely at the nominal energy of the 
resonance is apparent. Similar results are obtained for 
the 3/2 + and 3/2 _ resonances. 

According to the stabilization method, the eigenfunc- 
tions corresponding to the stabilized energies should cor- 
respond to well localized states, as expected for a res- 
onant state. This is confirmed in Fig. @] for the three 
resonances discussed above. In each panel, we compare 
the radial components of the scattering wave functions 
evaluated at the nominal energy of the resonance (solid 
lines), with the THO eigenfunction associated with the 
stabilized eigenvalue, using a basis of N=10 oscillator 
functions (dashed lines). Because the continuum wave- 
functions are not square-integrable, these functions have 
been conveniently scaled for a better comparison with 
the PS functions. 

Note that, for these three resonances, the channels cor- 
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FIG. 3: (Colour online) Eigenvalues obtained from the di- 
agonalization of the 11 Be Hamiltonian in a THO basis, as a 
function of the LST continuum parameter (7) in the upper 
panel, and as a function of the number of oscillator states 
included in the basis in the lower panel. The dashed line 
indicates the energy of the 5/2 + resonance and the dotted 
line, the energy of the 10 Be(2 + )+n threshold. 



responding to / = 2 are effectively bound, since the en- 
ergy of these resonances is below the n+ 10 Be(2 + ) thresh- 
old. The component based on the 10 Be(g.s.) is unbound 
but it shows the anticipated localization reminiscent of 
a quasi-stationary state. We note that, unlike the case 
of the bound states, we do not expect a perfect agree- 
ment between both calculations due to the exponential 
behaviour of the PS basis at large distances. Apart from 
that, it is also seen that, in the interior region, the four 
radial components are in very good agreement with the 
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FIG. 4: (Colour online) Radial parts of the continuum wave 
functions for the n Be resonances at e=1.2 MeV (5/2 + ), 2.7 
MeV (3/2") and 3.2 MeV (3/2+). The solid lines are obtained 
by direct integration of the Schrodinger equation, whereas 
the dashed lines are the result of the diagonalization in a 
THO basis with N=10 (5/2+) or N=9 (3/2", 3/2+) oscillator 
functions. For a better comparison, the normalization of the 
scattering state has been chosen in order to have the same 
magnitude as the discrete solution at the maximum. 



exact solution. 

The stabilization method provides also expressions for 
the width of the resonances in the PS basis [20|. How- 
ever, these expressions were originally developed for the 
single-channel case, and hence they cannot be directly 
applied to our case. To have an estimate of the width 
of the resonance we make use the density of states, de- 
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FIG. 5: Density of states for the 5/2+, 3/2+ and 3/2" reso- 
nances of the 11 Be nucleus, using different values of the basis 
size N. 



fined according to Eq. (fTTj) . This function is shown in 
Fig. E]for the discussed resonances, using different values 
of the basis size (N) . It can be seen how the density in- 
creases as more channels are open above the excitation 
energy of the core. It can be seen also that the pres- 
ence of a resonance gives rise to a peak in the density 
distribution. Based on this property, we have estimated 
the width of the resonance from the FWHM of the corre- 
sponding peak in the density distribution. For the 5/2 + , 
3/2~ and 3/2 + resonances considered above, this method 
yields T = 125 keV, 40 keV and 140 keV respectively. 
This widths are to be compared with the values reported 
in [11], namely, T = 125 keV, 50 keV and 100 keV. Except 
for the latter, for which our prescription gives a width 
40% larger, the agreement between both methods is very 
good in the other two cases. 

Just to complete our study, we show in Fig. |B] the 
comparison of the radial parts obtained by integration of 
the Schrodinger equation (solid lines) and by diagonal- 




FIG. 6: (Colour online) Radial parts of the 3/2~ wave func- 
tion for the n Be nucleus at E re i = 0.628 MeV obtained by 
direct integration of the Schrodinger equation (solid lines) and 
by diagonalization in a THO basis with N=15 states. All the 
components except the main one are multiplied by a factor of 
10. 



ization in a THO basis with N=15 (dotted lines) for a 
non-resonant state in the continuum. It can be observed 
that the agreement is also good for these states. 



B. Electric reduced transitions probabilities 

The electric transition probabilities provide also a use- 
ful test to assess the quality of the basis to represent the 
continuum states. These transition probabilities can be 
calculated using either the exact scattering states, using 
Eq. p^|) . or the pseudostates, using Eq. (J2T1) . In the lat- 
ter case, one obtains a discrete distribution, which can 
be converted to a continuous distribution by means of 
Eq. (f2"3"l) . In actual calculations, this equation is evalu- 
ated with a finite number of states (N) and hence this 
formula is only approximate. The degree of agreement 
of this approximate formula with the exact calculation 
provides a measurement of the quality of the PS basis 
to represent the continuum for a given operator. In this 
section we perform this test for the El and E2 operators. 

According to Eq. (TT6"|) . the electric operator for a va- 
lence+core system will contain in general contributions 
coming from the valence excitation, the core excitation 
and mixed excitations. However, in our test case, n Be, 
with core states restricted to the ground state (0 + ) and 
the first excited state (2 + ), dipole transitions will consist 
of pure single particle excitations. On the other hand, 
quadrupole transitions will contain both single particle 
and core excitations, but not simultaneous transitions. 
These simultaneous transitions will only affect octupole 
and higher order transitions, which will not be considered 
here. 

In Fig. [71 the energy distribution of the B(E\) obtained 
with a THO basis with N=20 functions is shown for 
n Be. Separate contributions for l/2~ and 3/2~ states 
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FIG. 7: (Colour online) Dipole electric transition probabil- 
ity (B(E1)) obtained with the THO basis and compared with 
experimental data from RIKEN [37| and MSU [jjg]. A THO 
basis with N=20 was used in the calculation, for which the 
calculated distribution is fully converged and indistinguish- 
able from the exact result using the exact scattering wave- 
functions. 



are shown by dotted and dashed lines, respectively. With 
this basis size, the calculated THO distributions are al- 
most indistinguishable from the exact calculation, ob- 
tained with the exact scattering states, so the latter has 
not been included in the figure. The available experi- 
mental distributions from two experiments performed at 
RIKEN Hzl and MSU J38| are also shown in the plot. 
The theoretical distribution lies in between the two ex- 
perimental sets of data. However, one has to keep in mind 
that the RIKEN data are inclusive with the respect to 
the 10 Be state and hence it might contain contributions 
where the core is left in an excited state. Moreover, it 
is also worth noting that the calculation will be sensitive 
to the choice of the n Be Hamiltonian. We have not ex- 
plored in this work this dependence since the purpose of 
this calculation is to test the quality of the basis, rather 
than a detailed comparison with the data. 

From Fig. [7] on sees that the calculated distribution 
shows a dip around e = 2.8 MeV, which is also visible 
in the data from Ref. [38[ . This behavior arises from the 
presence of the 3/2 _ resonance at this excitation energy. 
This resonance is relatively narrow (T = 50 keV) but 
it is only weakly coupled because it is mainly built on 
the excited core ( 10 Be(2 + )), whereas the ground state is 
mostly 10 Be(0+). 

Because only the single-particle excitation term of 
Eq. (|16j) contributes to this dipole transition, this observ- 
able can be also well reproduced within a single-particle 
model of Be, with the 10 Be core in its ground state, 
and including the appropriate spectroscopic factor for the 
| lo Be(0 + ) <g) 2s\/2) configuration. A departure from this 
behaviour is the afore-mentioned reduction of the B{E1) 
around 2.8 MeV which is due to a core dominated 3/2~ 
resonance. 

We have also evaluated the quadrupole electric transi- 
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FIG. 8: (Colour online) Quadrupole electric transition 
probability (B{E2)) obtained with the THO basis with 
N=20 oscillator functions. The upper panel is the calcu- 
lation including only the valence excitations. The bottom 
panel includes both valence and core contributions. 



tion probabilities, which are shown in Fig. [8] The dotted 
and dashed lines are the contributions from 3/2 + and 
5/2+ , respectively, whereas the solid line is the sum of 
both contributions. According to Eq. (fT6|) . in addition to 
the single-particle excitations, in this case we have also 
a contribution due to E2 transitions of the core which, 
in fact, give the main contribution to the total B(E2) 
strength. To illustrate better the contribution coming 
from the valence excitation and the core, we show in the 
upper panel of this figure the single-particle contribution, 
whereas in the bottom panel we show the full calcula- 
tion, including also contributions from the core. It is 
seen that the B(E2) strength is dominated by the core 
excitations, as expected for a collective transition. The 
peaks at e ~ 1.2 MeV and and e ~ 3.2 MeV are due to 
the 5/2+ and 3/2+ resonances. Unfortunately, no exper- 
imental or theoretical B(E2) for 11 Be has been found in 
the literature in order to compare with. 
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Application to the Coulomb breakup of 

208p b 



'Be on 



A more recent measurement of the Coulomb breakup 
of n Be can be found in the work by Fukuda and collab- 
orators [39|, who measured the breakup of a 11 Be beam 
at 69 MeV/nucleon on carbon and lead targets. 

At these energies and for very small angles the breakup 
is dominated by the Coulomb interaction. For angles 
below the grazing angle the differential break up cross 
section can be calculated semiclassically using the equiv- 
alent photon method (EPM) [|o[. For the El, which is 
expected to be the dominant one, the breakup cross sec- 
tion in the EPM method reads, 



\dndej bu 



16ir 3 dB(El)dN E1 (9 cm ,E x ) 



9hc de 



dQ r 



(24) 



where NEi(0 cm , E x ) denotes the number of virtual pho- 
tons with energy E x at scattering angle cm . In this 
treatment, the scattering angle corresponds to a classical 
Coulomb trajectory. The photon energy would be always 
E x = e + Sin- 

In a similar way, the E2 contribution to the breakup 
cross section which in this formalism is related to the 
B(E2) distribution, 



d 2 a 
dVtde 



4tt 3 
75hc 



Tic 



dB(E2) dN E2 
de dil cm ' 



(25) 



This contribution should be added to the dipole Coulomb 
break up. The equivalent photon number for E2 transi- 
tions can also be found in [401 ] . 

We have evaluated these contributions using the B(E1) 
and B(E2) distributions obtained with the THO ba- 
sis. Indeed, these expressions could be directly evalu- 
ated with the scattering states, since no discretization 
is required in this case. It is nevertheless illustrative to 
compare both calculations, to show the convergence of 
these observables with the size of the THO basis. In 
more sophisticated reaction models, such as the CDCC 
method with core excitation [3j, the use of a discretiza- 
tion method is mandatory, and hence the use of a discrete 
basis, like the THO proposed here, is more justified. 

From the expressions (|2"4"]l and (j2"5| , the energy and 
angular differential breakup cross sections are calculated 
by integrating in the scattering angle or in the excitation 
energy, respectively. In the former case, a critical ingredi- 
ent of the calculation is the minimum impact parameter 
(&min)- The model assumes that pure Coulomb breakup 
occurs only for b > 6 m ; n . By contrast, for b < b mm , the 
model assumes that other reaction mechanisms rather 
than pure Coulomb scattering take place (such as nuclear 
effects). Since these effects are not properly described by 
Eqs. (j2~4")) and (j2"5l) , these expressions are only evaluated 
for b > 6 m i n . This is also the reason why these Coulomb 
breakup experiments are focused at small angles, ideally 
below the grazing angle. In all the calculations the mini- 
mum impact parameter is settled according to the choice 
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FIG. 9: (Colour online) Energy (top panel) and angular 
(bottom panel) distributions of the exclusive breakup for 
the reaction u Be+ 208 Pb at 69 MeV/nucleon. The circles 
are the experimental data from Ref. [39(. The dot-dashed 
and dashed lines correspond to the El and E2 contribu- 
tions, calculated within the EPM method, and the solid 
line is the sum of both contributions. 



done in [39|, where the (El) Coulomb breakup cross sec- 
tion is also evaluated using the EPM method. 

In Fig. [9] we compare the calculated energy (upper two- 
panel figure) and angular differential cross sections (lower 
panel) with the experimental data. The separate El and 
E2 contributions, as well as their sum, are shown in each 
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panel. The calculations have been convoluted with the 
ex per imental angular and energy distributions reported 
in [33] . It is clearly seen that the main contribution comes 
from the dipole break up. In the angular distribution, 
the sum of both contributions cannot be distinguished at 
the smaller angles from the pure El contribution. The 
small E2 contribution is only observed in the energy re- 
gions of the resonances. This difference is nevertheless 
washed out once the energy resolution of the experiment 
is considered. Despite this small contribution, it is ob- 
served that the E2 component improves the agreement 
in the energy region nearby the 5/2 + resonance. Com- 
paring the total and dipole angular distributions of the 
Coulomb break up one can infer up to what angle one 
should consider pure El excitations in order to extract a 
B(E\) distribution not affected by the 5/2 + resonance. 



V. SUMMARY AND CONCLUSIONS 

We have investigated the problem of the description of 
the states of a particle moving in a deformed potential in 
terms of a pseudo-state (PS) basis. In the PS method, the 
states of the system are approximated by the eigenstates 
of the Hamiltonian in a basis of square-integrable func- 
tions. The negative eigenvalues arc identified with the 
bound states of the system, whereas the positive eigen- 
values are regarded as a discrete and finite representation 
of the continuum spectrum. Identification of resonances 
is done using the so-called stabilization method [13, [20| • 

Following our previous choice for non-deformed sys- 
tems, we propose to use as PS basis the Transformed 
Harmonic Oscillator (THO) basis. The basis functions 
are obtained by applying an analytic local scale transfor- 
mation [3, U3, to the conventional HO basis. The trans- 
formation is such that it converts the Gaussian asymp- 
totic behavior of the HO function into an exponential. 

The method has been applied to the n Be nucleus, 
treated within a particle-rotor model. The 10 Be core 
is assumed to have a permanent axial deformation with 
P% = 0.67 [lg]. We have shown that the bound-state en- 
ergies and wave functions are very well described using 
a relatively small basis, showing perfect agreement with 
those obtained by direct integration of the Schrodinger 
equation. We have shown that the resonances 5/2 + , 3/2~ 
and 3/2 + are also well described with the method using 
small THO bases. It has also been checked that the wave 
functions of the non-resonant continuum calculated with 
the THO method compare well with the state computed 



by direct integration of the Schrodinger equation at the 
same energy. 

We have given expressions for the El and E2 electric 
transition probabilities in the discrete basis, and we have 
proposed a method to obtain smooth distributions from 
these discrete values. To illustrate this method, we have 
calculated the B(E1) and B(E2) electric transition prob- 
abilities for the n Be nucleus. These distributions show 
a fast convergence rate with the basis size, and the con- 
verged results are in perfect agreement with the exact cal- 
culation, obtained with the exact scattering states. With 
the adopted Hamiltonian, the calculated B(E1) distribu- 
tion is consistent, but somewhat larger, than the experi- 
mental data from MSU 38] . 

Finally, we have applied the model to the Coulomb 
breakup of n Be on 208 Pb at 69 MeV/nucleon, comparing 
with the data from Ref. [33]. The reaction is treated 
in a semi-classical picture, using the equivalent photon 
method, and including both El and E2 contributions. 
The calculated angular distribution is in good agreement 
with the data for scattering angles below 3°. Beyond this 
angle, other effects not considered in the EPM method, 
such as nuclear breakup, are expected to take place. The 
calculated energy distribution is also in good agreement 
with the data, particularly when the angular range is 
below the grazing angle. 

We conclude that the THO basis provides a suitable 
representation to describe two-body composite systems 
(bound and unbound states) including the core defor- 
mation. This study provides the needed test for accom- 
plishing a similar study for more interesting cases, such 
as three-body composite systems including core deforma- 
tion or three-body scattering problems (two-body pro- 
jectile plus a target) including dynamic core excitation. 
Work toward this direction is in progress. 
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